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Abstract. We have calculated the density-density (Lindhard) response function of 
a homogeneous two-dimensional (2D) hole gas in the static (w = 0) limit. The 
' ^ I bulk valence-band structure comprising heavy-hole (HH) and light-hole (LH) states is 

modeled using Luttinger's k-p approach within the axial approximation. We elucidate 
CJ i how, in contrast to the case of conduction electrons, the Lindhard function of 2D holes 

exhibits unique features associated with (i) the confinement-induced HH-LH energy 
^S| ■ splitting and (ii) the HH-LH mixing arising from the charge carriers' in-plane motion. 
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1. Introduction 

The density-density response function is a very fundamental materials characteristics, 
as it determines a host of thermodynamic and transport properties in condensed- 
matter systems [1]. It has been discussed extensively within the paradigmatic model of 
the homogeneous electron gas [2] and studied for low-dimensional conductors realised 
in semiconductor heterostructures [3]. More recently, static and dynamic response 
properties of two-dimensional (2D) conduction-electron systems with spin-orbit coupling 
have been investigated in considerable detail [4-9]. This surge of interest arose partly 
because of important ramifications for possible spintronics applications [10,11]. In 
contrast, very few studies have considered how the peculiar electronic properties of a 
typical semiconductor's valence band [12] affect the polarizability and other many-body 
response functions of p-tjpe semiconductor materials, and these existing works [13-15] 
have focused on bulk (3D) systems. As high-quality 2D holes gases have recently become 
available for experimental study, both in modulation-doped [16-20] and accumulation- 
layer [21] heterostructures, a detailed theoretical analysis of their many-body response 
properties is warranted. Here we provide such a study and show how the intricate 
interplay between quantum confinement and strong spin-orbit-coupled dynamics in the 
valence band [22] has a profound effect on the static polarizability. 

Charge carriers from the conduction and valence bands of typical semiconductors 
exhibit profoundly different spin properties. Conduction electrons are quite ordinary in 
that they are spin- 1/2 particles carrying a fixed intrinsic magnetic dipole moment, like 
free electrons in vacuum. Holes are different; they have an intrinsic spin-3/2 degree of 
freedom because valence-band states are strongly modified by spin-orbit coupling [12]. 
As a result, the orbital dynamics of holes in a bulk sample also depends on the magnitude 
of projection for their spin parallel to their direction of motion. States with spin-3/2 
projection quantum number mj = ±3/2 (±1/2) are called heavy holes, HHs (light holes, 
LHs), because their band-energy dispersion has a smaller (larger) curvature. When holes 
are confined in a 2D heterostructure, the quantum-well growth direction is the natural 
spin-quantisation axis (taken to be the z direction in the following), and the difference 
in effective masses translates into an energy splitting between the HH and LH subband 
edges corresponding to the same transverse orbital bound state [23]. As the in-plane 
motion couples HH and LH states, 2D holes with finite wave vector ky = {kx,ky) are 
no longer of purely HH or LH type [22,23]. While the HH-LH (subband-edge) splitting 
is easily accounted for and usually included in theoretical analyses, the HH-LH mixing 
has sometimes been ignored. It may be tempting to make such a simplification, given 
that the density of typical 2D hole gases is quite often low enough that only the lowest 
(HH-like) subband is occupied. However, detailed analysis shows this approach to be 
too crude for most relevant situations [22]. Even only qualitatively accurate predictions 
basically always require the inclusion of HH-LH mixing alongside the HH-LH splitting. 
As we will show below, the density-density response of 2D hole gases is strongly affected 
by HH-LH mixing, i.e., is not simply the sum of the response functions of independent 
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2D (HH and LH) gases. 

This article is organised as follows. In Section 2, we introduce our model for the 
upper-most valence band of typical semiconductors, which is based on the Luttinger 
Hamiltonian [24,25] in axial approximation [22,26,27]. The definition and basic 
calculational details for the density-density (Lindhard) response function are given in 
Sec. 3, including analytical results pertaining to the 2D hole gas in certain limits. We 
present plots of the numerically determined static polarizability in Sec. 4 and discuss 
basic features. Our conclusions are given in Sec. 5. 

2. Luttinger-model description of a 2D hole system 

The Luttinger model [25] provides a useful description of the upper-most valence band 
of typical semiconductors in situations where its couphngs to the conduction band and 
split-off valence band are irrelevant. We adopt this model here to investigate how 
the many-body physics of 2D holes is affected by their peculiar spin-3/2 properties. In 
principle, more extended [27, 28] multiband Hamiltonians could be employed to improve 
the accuracy of quantitative predictions. However, to illustrate the qualitatively new 
features exhibited by 2D hole gases in contrast to their conduction-electron counterparts, 
the Luttinger-model description is adequate. The particular geometry for our case of 
interest suggests using, as our starting point, the Luttinger-model Hamiltonian Hi^ in 
axial approximation [22,26,27]: 

Hi^ = Ho + H, + H2 , (la) 

H,= -^ ^,iki + I,i)+^,i\^i-2k-i]iJ^--l] , (lb) 



2mQ 



7i(kJ + A:,^)+7i(k5-2fc9(jf-^i) 



mo \ / 

Cartesian components of the spin-3/2 matrix vector are denoted by Jx,y,z, and we used 
the abbreviations k± = k^ i: iky, J± = [J^ ± iJy)/\/2, and {A,B} = (AB + BA)/2. 
The constants 71 and 7j are materials-dependent bandstructure parameters [29]. Note 
that the 7^ depend also on the quantum-well growth direction; their explicit expressions 
in terms of the standard Luttinger parameters [25,29] 72 and 73 can be found, e.g., in 
Table C.lOof Ref. [22]. 

The dynamics of holes confined in a 2D quantum well is modeled by the Hamiltonian 
Hi^ + V{z). In the following, we assume the external potential V{z) to be a hard- 
wall confinement of width d and consider only its lowest size-quantised orbital bound 
state. An effective Hamiltonian describing the 2D hole gas is then obtained from (la) 
by replacing kz — )■ (kz) = and /c^ — )■ (A;^) = (tt/c?)^. Introducing the energy scale 
Eq = 7r^/i^7i/(2morf^) and measuring wave vectors in units of n/d, the 2D hole-gas 
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Hamiltonian is given by 



H, 
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Here k^.y = kx^yd/n, 7 = 71/71, and a = 73/71. We are using the parameterisation in 
terms of 7 and a to be able to separately discuss the effects of HH-LH splitting, which 
is embodied in Hq , and HH-LH mixing arising from i/mix- 

Diagonalising H^ from (2a) yields in-plane dispersion relations £'j(k||) = —E^ e^ 
with j = 1, 



4 and 



Sj) 



l + k; + k 



Here (Ti = o"2 = — o"3 
Fermi wave vectors 

kp.o = 



y + ^.7 V (^^ + ^^ - 2)^ + 3a2(A;2 + A;2)2 . (3) 

o"4 = 1. Using the result (3), we obtain the two dimensionless 



£f - 1 - 27' T 7V(^F-3)2 + 3a2[(ep-i)2_4^2] 
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in terms of the dimensionless Fermi energy ep = —Ey/Eq. The 2D-hole sheet density 
n2D is related to the dimensionless Fermi wave vectors according to 

-Try kl^Qjsp - [1 + 27]) + kl^Qjep - [1 - 27]) 
.d) 2tx 

where Q{x) denotes the Heaviside step function. 

The eigenvectors |Xk ) corresponding to eigenvalues ii^j(k||) of H\^ can be 

straightforwardly determined. For j = 1,2 and in the basis representation where Jz 
is diagonal, we find 



(i)\ 
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(3)\ 



(6) 
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with s = kji — 2 and t = \ s"^ + Sa^kn. The remaining eigenspinors 1x^7 ^"^^ IXk 7 ^^^ 



71) 



(2)\ 



obtained by changing t — )■ — t in |Xk ) S'-^^d |Xk )) respectively. As the scalar products 
(Xk IXk +q) enter in the calculation of the Lindhard function, we briefly discuss their 
relevant properties. The moduli |(Xk IXk +q)l ^^^ found to be equal for all j = 1, . . . ,4. 
Also, |(Xk IXk +q)l ^^e pairwise the same for (j, /) = (1,3) and (3, 1); and (j, /) = (2,4) 
and (4, 2). These relations can be verified using the explicit form of the column vectors 
in (6) together with the fact that the eigenspinors satisfy orthonormality relations. 
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3. Lindhard function of a 2D hole gas: General expression and special cases 

The general definition [2] of the Lindhard function, specialised to a 2D system, reads 

^(u a) - lim V f "^'^'1 Kv^-^'^lvW M2 np[E,{k\\)] - n^jEim + q)] 
X(u;,q)-hm2^y (^K^k„ IXk„+c)l /^ + ,5 + i?,.(k,|) - i?,(k,| + q) ' ^^ 

with npiE) denoting the Fermi-Dirac distribution function. The expression given in (7) 
can be simplified by using a polar-coordinate representation where k^ = k\\ cos0 and 
ky = /cii sin and performing a change of variables in the terms involving nF[Ej{\i.\\ + q)] 
such that k|| — )■ ky — q and — )■ + tt, which leaves the energy difference and the spinor 
overlap invariant. Using also the description in terms of dimensionless quantities and 
specialising to the zero-temperature limit, the Lindhard function can be expressed as 
X(w,q) = -(2mo//i^7i)x(^,g), with 






X{uj, 9) = lim 5^ 5Z / ' ^11 ^^11 / 

», — LI .' i_i "^ Jo 



We use a notation where Kp^ = Kp^ = kp^, Kp.^ = Kp^ = kp^, and u = hu/Eo. 
Note that, because the Luttinger Hamiltonian in axial approximation exhibits rotational 
invariance of in-plane hole motion, the Lindhard function depends on wave vector q only 
via its (dimensionless) magnitude q. Also, within our effective 2D description, xi^^ q) 
is independent of the quantum- well width d and inversely proportional to 71. 

In the following, we consider the static limit, which is obtained by setting a) = 0. 
Specialising further to certain limiting situations, we can find analytical expressions for 
the Lindhard function. For example, for the case of vanishing HH-LH mixing obtained 
by letting a — > 0, the matrix |(Xk IxL +q)P ^f modulus-squared scalar products reduces 
to the unity matrix, and the simple analytical expression 



is found, where 0"i = — (T2 = 1- Inspection of the result (9) shows that, with only HH- 
LH splitting included, the static Lindhard function comprises two separate HH and LH 
contributions, each being the standard 2D-electron-gas expression [2] with Fermi wave 
vector and effective mass adjusted to the respective HH and LH values. On the other 
hand, taking the limit q — ?■ in (8), the matrix of modulus- squared spinor overlaps 
again becomes the unity matrix, and we find an analytical result for the (dimensionless) 
density of states at the Fermi energy. 



.7=1 



1 - a,-7 ■ 



(2-A;2,)2 + 3a2A;4 



(10) 



Thus we see that one effect of HH-LH mixing is to introduce an energy (and concomitant 
density) dependence into the density of states of 2D holes. 
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4. Static polarizability of 2D holes: Numerical method and results 

With analytical expressions unavailable for the Lindhard function (8) in the more general 
case with both q and a finite, we have to resort to numerical calculations to investigate 
in greater detail how HH-LH mixing affects the static polarizability x(0)^)- Note that 
the latter is an entirely real-valued function. The procedure for its numerical calculation 
is explained in the following Subsection, and our results are given thereafter. 

4.I. Brief outline of the calculational method 

For u = 0, the integrand of (8) has poles whenever the energy difference in the 
denominator vanishes. These poles are regularised by the parameter S, which needs 
to be set to zero after performing the integrations. We calculate these integrals 
numerically, taking special care in the regions close to values of the integration variables 
corresponding to a vanishing denominator. To identify the pole structure of the 
Lindhard function, we write the inverse of the energy difference as 

+ ((^i.S + S,,4)i6i,3 + 5l,4)^^f^, (11) 

where 6j^i denotes Kronecker's delta symbol. The quantities appearing in (11) are 

ai,2 = Tl\/4: - AB + (1 + 3a^)B + f + k\\qcos(t) (12a) 
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6 = 7W(-2 + A;2 + g2 + 2A;||gcos0)2 + 3«2(A;J + f + 2k\\qcos(py . (126) 
The denominators in (11) can be written as 

(13) 



a? 2 - 62 4^2g2[l _ (1 + 3^,2)^2] X^^^ - Y,^2 \COS (f) - Xi,2 COS - 11,2 

with the positions of the poles given by 



Xi: 



2h 



±27./4 - 4A;2 + (1 + 3a2)fc4 _ ^2 ^ ^2^ ^ g^2)^2 _ 4 + ^2 ^ 3^2^2] 

Yi 2 = ^ . (14) 

' 2fc||g[l-(l + 3a2)^2] ^ > 

As can be seen from (13), poles are encountered in the integration over when 
|Ari,2|, |li,2| < 1- We have employed a Cauchy principle-value integration to regularise 
the Lindhard function in the vicinity of the poles specified in (14). 
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Figure 1. (a) Normalised Lindliard function xiQ) for tire static limit and (b) the 
quantity Ax{q,a) that measures the impact of HH-LH mixing (see text), plotted as 
a function of wave-vector magnitude q. The blue curve is for 7 = 0.31 and a = 1.2, 
which are parameter values applying to an [001]-grown heterostructure in GaAs. For 
comparison, results are also shown for a = 1 (red curve), 0.5 (magenta curve), and 
(green curve). For all cases, a value n = 0.0608 for the dimensionless 2D hole density 
was used. 



4-2. Results for model parameters applying to an [001] quantum well in GaAs 

High-quality 2D hole gases have recently been fabricated from [001]-grown GaAs 
heterostructures [16,17,20]. To obtain results applicable to these systems, we use the 
appropriate model parameters 7 = 0.31 and a = 1.2. Results for this configuration are 
presented below. For comparison and to clearly show the impact of HH-LH mixing, we 
also show results for the static polarizability when a = 1, 0.5, and 0. In all these cases, 
we limit ourselves to the low-density regime where only the highest, HH-like, 2D subband 
is occupied. The reason for this caution is the fact that, within our model using a hard- 
wall confinement, the spectrum of all other than the highest subband poorly matches 
that of the real GaAs sample. To be specific, we choose n = 0.0608. Recalling that the 
2D-hole sheet density is related to the dimensionless density by n = {-k /dYn, this value 
would correspond to a density of n = 1.5 x 10^^ m~^ in a 20-nm quantum well. 

To avoid cluttering our notation, we suppress the zero-frequency argument in the 
formal expression of the static Lindhard function from now on: ^(O^q) — xio)- Iii 
Figure 1(a), we plot x(Q')/x(0) as a function of q/y/STm, for different values of the 
parameter a that quantifies the HH-LH mixing. It is apparent that a finite a leads to a 
significant suppression of xil) below the constant-plateau value usually associated with 
2D systems [2] for q < VSvm. To make the impact of HH-LH mixing quantitatively 
explicit, we define the variable 



Ax(g;a) 



1- 



x{(l) 



X(0) 



xiQ) 



X(0) 



-1 



(15) 



where xIq) is the analytical result (9) obtained for the limit a = but with Fermi 
wave vectors adjusted to coincide with those found in the case of the finite a under 
consideration. Thus the function Ax(g; a) measures the relative change exhibited in 
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Figure 2. Normalised Lindliard function x(g) for tlie low-density regime of a fictitious 
semiconductor material with 7 — 0.2 and (dimensionless) density n — 0.0608. The red, 
magenta, and green curves correspond to a = 1, 0.5, and 0, respectively. 



the normalised static polarizability that is due to a finite a but goes beyond a simple 
renormalisation of Fermi wave vectors. | In Figure 1(b), we show Ax(g; a) for a = 1 (red 
curve) and a = 1.2 (blue curve). It shows a strong variation as a function of g/yovm 
and reaches the 10% level. 

4-3. Results for a model semiconductor: High- and low-density regimes 

For high-enough 2D densities, holes will occupy both the HH-like and LH-like subbands 
arising from the lowest-energy orbital bound state in the quantum well. It can be 
expected that this high-density regime is qualitatively different from the situation at 
low density where only the highest (HH-like) 2D subband is occupied. To treat the 
case of high density consistently within our adopted model, it needs to be ensured that 
the LH-like subband arising from the lowest-energy orbital bound state is still higher in 
energy than the HH-like subband associated with the next orbital-bound-state level. For 
a hard-wall confinement considered here, a system with 7 = 0.2 satisfies that condition. 
Although this value does not directly correspond to a specific semiconductor material, 
we use it to illustrate the generically different impact of HH-LH mixing in the low and 
high-density regimes, respectively. 

To provide a clear benchmark for comparing high and low-density regimes, we start 
by presenting the result for the low density case (n = 0.0608, same value as used in the 
calculations for Figure 1) in Figure 2. The obtained curves look qualitatively similar 
to those found for the low-density regime in GaAs (different 7, shown in Figure 1), but 
the quantitative level of suppression below the plateau value obtained in the hmit of 
vanishing a is different here. 

I In the low-density limit considered here, there is only one Fermi wave vector whose magnitude is the 
same for all values of a, and xil) actually coincides with x{l)\a=o- However, as we will see further 
below, xil) 7^ x('?)|q=o in the more general case when both the HH and LH subbands arising from the 
lowest 2D orbital bound state are occupied and, thus, two Fermi wave vectors exist. 
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Figure 3. (a) Normalised Lindhard function x{q) for the static limit and (b) the 
quantity Ax{q',a) that measures the impact of HH-LH mixing (see text), for the case 
with 7 = 0.2 and (dimensionless) density n = 0.4055, corresponding to the high-density 
regime with both HH-like and LH-like 2D subbands occupied. The red, magenta, and 
blue curves are obtained for a = 1, 0.5, and 0, respectively. For comparison, we 
also plot (as the green curve) the analytical result (9) with Fermi-wave-vector values 
adjusted to coincide with the a = 1 case. Therefore, the red and green curves in 
panel (a) exhibit kink-like features at the same values of q, and deviations between 
the two illustrate the effect of HH-LH mixing beyond a simple renormalisation of the 
Fermi wave vectors. 



As an illustration of the high-density regime, we present results for n = 0.4055, 
which would correspond to n = 10^^ m~^ in a 20-nni quantum well. The normalised 
Lindhard function for this case is plotted in Figure 3(a). Note that, for the three values 
of a for which results are presented, the two Fermi wave vectors are different, see (4). 
As a result, the sharp features arising in the Lindhard function from poles at —q/2kp^ ^ 
appear at different values of q in each curve. In contrast to the low-density case, a 
plateau is again exhibited in the static Lindhard function (for q < 2A;fi)- To illustrate 
effects due to HH-LH mixing beyond a simple renormalisation of the two Fermi wave 
vectors, we also show (as the green curve) the analytical result (9) for the Lindhard 
function in the limit a = but with values for the Fermi wave vectors taken from the 
case a = 1. The latter result corresponds to that expected for two independent 2D hole 
gases with different Fermi wave vectors. Thus, any deviation between the red and green 
curves is entirely due to the mixed HH-LH character of 2D hole states. Figure 3(b) 
shows the corresponding difference function Ax(g; a) for a = 1. Again, HH-LH mixing 
appears to cause differences around the 10%-level. 



5. Summary and Conclusions 

We have calculated the density- density response function of a homogeneous 2D hole 
gas in the static hmit, based on the Luttinger-model description of the upper-most 
valence band within the axial approximation. While this approach neglects the warping 
of energy dispersions (and, hence, Fermi surfaces) for the holes' in-plane motion due to 
the cubic crystal symmetry, it already captures essential new features arising from the 
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peculiar valence-band properties. We furthermore focused only on the lowest orbital 
bound state in a symmetric quantum well defined by a hard- wall confinement. HH-LH 
splitting gives rise to the existence of two energetically separated 2D hole subbands, 
one (at higher energy) mostly HH-like and the other of mostly LH character. However, 
except for states with zero in-plane kinetic energy, HH and LH amplitudes are mixed, 
and we have elucidated how this mixing gives rise to marked changes in the shape and 
magnitude of the static density response function. New analytical results are derived 
for the limit g — )■ 0, but the case with finite q and HH-LH mixing included could only 
be treated numerically. 

Both for practical and theoretical reasons, it makes sense to distinguish two basic 
situations. One corresponds to the low-density regime where only the highest (mostly 
HH-like) 2D subband is occupied. In this limit, our model can be expected to describe 
real semiconductor heterostructures quite accurately, even quantitatively, if adequate 
band-structure parameters are used as input to our calculations. As it turns out, even 
though holes are present only in the HH-like 2D subband, HH-LH mixing importantly 
affects the static density response. In particular, the response function is suppressed 
below the plateau exhibited by the standard 2D-electron-gas result [2], as illustrated 
in Figures 1(a) and 2. In the high-density regime, where both the HH-like and LH- 
like 2D subbands associated with the lowest-energy quantum-well bound state are 
occupied, the density response differs from that expected for two independent 2D 
hole gases. Thus HH-LH mixing is shown to influence the Lindhard function beyond 
a trivial renormalisation of Fermi-wave-vector magnitudes. We have defined, and 
calculated, the quantity Ax(g; a) to make the nontrivial effects arising from HH-LH 
mixing quantitatively explicit. For the parameters considered, relative changes on the 
order of 10% are seen. 

In this work, we employed the approximation of hard-wall (infinite-height) 
quantum-well barriers. In the more realistic case of finite-height barriers, the HH and 
LH bound-state wave functions will both penetrate into the barrier regions. Due to 
HH-LH splitting, the range of penetration will be higher for LH states, i.e., the effective 
quantum-well width will be larger for LHs. This subtle difference between HH and LH 
bound states scales with the parameter 7 and can be expected to result in corrections 
of order 07^ to effects due to HH-LH mixing. 

The qualitative and quantitative impact of HH-LH mixing can be expected to 
affect physical properties of 2D hole gases in an important way and, thus, render 
their behaviour quite different from that exhibited by 2D conduction-electron systems. 
Examples for physical observables that are affected include the shape and range of 
Friedel oscillations exhibited by 2D hole gases in response to impurity charges present, 
e.g., in the doping layer of modulation-doped heterostructures. HH-LH mixing should 
then also influence the 2D-hole-mediated Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction between magnetic impurities. A detailed investigation of this effect could 
shed new light on how to tailor the ferromagnetic properties of 2D diluted-magnetic- 
semiconductor heterostructures [30-34]. 
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